Carbonates in space - The challenge of low temperature data

Carbonates have repeatedly been discussed as possible carriers of stardust emission bands. However, the band assignments proposed so far were mainly based on room temperature powder transmission spectra of the respective minerals. Since very cold calcite grains have been claimed to be present in protostars and in Planetary Nebulae such as NGC 6302, the changes of their dielectric functions at low temperatures are relevant from an astronomical point of view. We have derived the IR optical constants of calcite and dolomite from reflectance spectra - measured at 300, 200, 100 and 10K - and calculated small particle spectra for different grain shapes, with the following results: i) The absorption efficiency factors both of calcite and dolomite are extremely dependent on the particle shapes. This is due to the high peak values of the optical constants of CaCO3 and CaMg[CO3]2. ii) The far infrared properties of calcite and dolomite depend also very significantly on the temperature. Below 200K, a pronounced sharpening and increase in the band strengths of the FIR resonances occurs. iii) In view of the intrinsic strength and sharpening of the 44 mum band of calcite at 200-100K, the absence of this band -- inferred from Infrared Space Observatory data -- in PNe requires dust temperatures below 45K. iv) Calcite grains at such low temperatures can account for the '92' mum band, while our data rule out dolomite as the carrier of the 60-65 mum band. The optical constants here presented are publicly available in the electronic database http://www.astro.uni-jena.de/Laboratory/OCDBComment: 20 pages, 10 figures, accepted by ApJ, corrected typo

From interstellar abundances to grain composition: the major dust constituents Mg, Si and Fe

We analyse observational correlations for three elements entering into the composition of interstellar silicate and oxide grains. Using current solar abundances (Asplund et al. 2009), we convert the gas-phase abundances into dust-phase abundances for 196 sightlines. We deduce a sharp difference in abundances for sightlines located at low (|b|<30\degr) and high (|b|>30\degr) galactic latitudes. For high-latitude stars the ratios Mg/Si and Fe/Si in dust are close to 1.5. For disk stars they are reduced to ${\rm Mg/Si} \sim 1.2$ and ${\rm Fe/Si} \sim 1.05$. The derived numbers indicate that 1) the dust grains cannot be the mixture of silicates with olivine and pyroxene composition only and some amount of magnesium or iron (or both) should be in another population and 2) the destruction of Mg-rich grains in the warm medium is more effective than of Fe-rich grains. We reveal a decrease of dust-phase abundances and correspondingly an increase of gas-phase abundances with distance $D$ for stars with D\ga 400\,pc. We attribute this fact to an observational selection effect: a systematic trend toward smaller observed hydrogen column density for distant stars. We find differences in abundances for disk stars with low (E({\rm B-V}) \la 0.2) and high (E({\rm B-V}) \ga 0.2) reddenings which reflect the distinction between the sightlines passing through diffuse and translucent interstellar clouds. For Scorpius-Ophiuchus we detect an uniform increase of dust-phase abundances of Mg and Si with an increase of the ratio of total to selective extinction $R_{\rm V}$ and a decrease of the strength of the far-UV extinction. This is the first evidence for a growth of Mg-Si grains due to accretion in the interstellar medium.Comment: 16 pages, 16 figures, accepted for publication in Astronomy and Astrophysic

Migration and Accretion of Protoplanets in 2D and 3D Global Hydrodynamical Simulations

Planet evolution is tightly connected to the dynamics of both distant and close disk material. Hence, an appropriate description of disk-planet interaction requires global and high resolution computations, which we accomplish by applying a Nested-Grid method. Through simulations in two and three dimensions, we investigate how migration and accretion are affected by long and short range interactions. For small mass objects, 3D models provide longer growth and migration time scales than 2D ones do, whereas time lengths are comparable for large mass planets.Comment: 4 pages, 4 figures; to appear in the Conference Proceedings of "Scientific Frontiers in Research on Extrasolar Planets

Do we need to know the temperature in prestellar cores?

Molecular line observations of starless (prestellar) cores combined with a chemical evolution modeling and radiative transfer calculations are a powerful tool to study the earliest stages of star formation. However, conclusions drawn from such a modeling may noticeably depend on the assumed thermal structure of the cores. The assumption of isothermality, which may work well in chemo-dynamical studies, becomes a critical factor in molecular line formation simulations. We argue that even small temperature variations, which are likely to exist in starless cores, can have a non-negligible effect on the interpretation of molecular line data and derived core properties. In particular, ``chemically pristine'' isothermal cores (low depletion) can have centrally peaked C$^{18}$O and C$^{34}$S radial intensity profiles, while having ring-like intensity distributions in models with a colder center and/or warmer envelope assuming the same underlying chemical structure. Therefore, derived molecular abundances based on oversimplified thermal models may lead to a mis-interpretation of the line data.Comment: ApJL, accepte

Fragmentation, infall, and outflow around the showcase massive protostar NGC7538 IRS1 at 500 AU resolution

Aims: Revealing the fragmentation, infall, and outflow processes in the immediate environment around massive young stellar objects is crucial for understanding the formation of the most massive stars. Methods: With this goal in mind we present the so far highest spatial-resolution thermal submm line and continuum observations toward the young high-mass protostar NGC7538 IRS1. Using the Plateau de Bure Interferometer in its most extended configuration at 843mum wavelength, we achieved a spatial resolution of 0.2"x0.17", corresponding to ~500AU at a distance of 2.7\,kpc. Results: For the first time, we have observed the fragmentation of the dense inner core of this region with at least three subsources within the inner 3000 AU. The outflow exhibits blue- and red-shifted emission on both sides of the central source indicating that the current orientation has to be close to the line-of-sight, which differs from other recent models. We observe rotational signatures in northeast-southwest direction; however, even on scales of 500 AU, we do not identify any Keplerian rotation signatures. This implies that during the early evolutionary stages any stable Keplerian inner disk has to be very small (<=500 AU). The high-energy line HCN(4-3)v2=1 (E_u/k=1050K) is detected over an extent of approximately 3000 AU. In addition to this, the detection of red-shifted absorption from this line toward the central dust continuum peak position allows us to estimate infall rates of ~1.8x10^(-3)Msun/yr on the smallest spatial scales. Although all that gas will not necessarily be accreted onto the central protostar, nevertheless, such inner core infall rates are among the best proxies of the actual accretion rates one can derive during the early embedded star formation phase. These data are consistent with collapse simulations and the observed high multiplicity of massive stars.Comment: Accepted for Astronomy & Astrophysics, 8 pages, also available at http://www.mpia.de/homes/beuther/papers.htm

Effective medium theories for irregular fluffy structures: aggregation of small particles

We study the extinction efficiencies as well as scattering properties of particles of different porosity. Calculations are performed for porous pseudospheres with small size (Rayleigh) inclusions using the discrete dipole approximation. Five refractive indices of materials covering the range from $1.20+0.00i$ to $1.75+0.58i$ were selected. They correspond to biological particles, dirty ice, silicate, amorphous carbon and soot in the visual part of spectrum. We attempt to describe the optical properties of such particles using Lorenz-Mie theory and a refractive index found from some effective medium theory (EMT) assuming the particle is homogeneous. We refer to this as the effective model. It is found that the deviations are minimal when utilizing the EMT based on the Bruggeman mixing rule. Usually the deviations in extinction factor do not exceed $\sim 5$ for particle porosity ${\cal P}=0 - 0.9$ and size parameters x_{\rm porous} = 2 \pi r_{\rm s, porous}/\lambda \la 25. The deviations are larger for scattering and absorption efficiencies and smaller for particle albedo and asymmetry parameter. Our calculations made for spheroids confirm these conclusions. Preliminary consideration shows that the effective model represents the intensity and polarization of radiation scattered by fluffy aggregates quite well. Thus, the effective models of spherical and non-spherical particles can be used to significantly simplify computations of the optical properties of aggregates containing only Rayleigh inclusions.Comment: 24 pages, 9 figures, accepted for publication in Applied Optic